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Abstract: 

Domain walls in equilibrium phase transitions propagate in a preferred direction so 
as to minimize the free energy of the system. As a result, initial spatio-temporal patterns 
ultimately decay toward uniform states. The absence of a variational principle far from 
equilibrium allows the coexistence of domain walls propagating in any direction. As a con- 
sequence, persistent patterns may emerge. We study this mechanism of pattern formation 
using a non-variational extension of Landau's model for second order phase transitions. 
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Second order phase transitions, such as in ferromagnets or liquid-vapor systems, are 
manifestations of spontaneous symmetry breaking occurring near thermal equilibrium. 
The coexistence of broken symmetry states beyond the transition point gives rise to spa- 
tial patterns consisting of domain walls or fronts separating regions of different phase. The 
dynamics of domain walls near thermal equilibrium, are dictated by a variational princi- 
ple, namely, the minimization of the free energy. As a consequence, initial spatio-temporal 
patterns ultimately decay toward the stationary homogeneous states of lowest free energy. 
Front structures are commonly observed in far-from-equilibrium systems as well. Walls 
separating conductive and convective states in binary mixtures [1], excited and recovery 
regions in autocatalytic chemical reactions [2], or different phase-locked states in paramet- 
rically forced surface waves [3], are a few examples. Unlike equilibrium systems, however, 
no general variation principle, like the minimization of free energy, applies for systems 
maintained far from equilibrium. The possible outcome, as emphasized recently [4-7], 
is the appearance of localized structures and, more generally, persistent spatio-temporal 
patterns. 

In this paper we further elaborate on pattern formation as a no n- variational effect. 
We consider spatio-temporal patterns involving domain walls and show that multiplicity 
of stable front solutions may give rise to persistent patterns. As a model system we choose 
to study a non-variational extension of the Landau- Ginzburg model for a scalar order 
parameter. The extended system takes the form of coupled reaction diffusion equations 
that have extensively been studied in the context of chemical and biological patterns [8-13]. 
Some of the results to be described here have already been obtained in that context before, 
particularly by Rinzel and Terman [12]. We rederive these results and present them here 
in a way that best illustrates the different point of view we put forward in this paper. 

Consider a variational system whose free energy (Liapunov functional) is given by 



where U{4>, h) = —(fi 2 /2 + 4> 4 /4: + h(fi, 4>{x,t) is a scalar order parameter, h is a constant 




(1) 
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bias field and the subscript x denotes the spatial partial derivative. For h values in the 
range -2/(3^) < h< 2/(3y/3) the free energy density, U, has a double-well form. The 
two wells correspond to stationary homogeneous states, characterized by order-parameter 
values (j)-(h) and 4>+(h) that solve the cubic equation (fi s — </> + h = 0. The relaxation 
toward any of these stationary states is governed by the equation 4> t = —SJ 7 / 8(f), or 

<h = f(<l>,h) + (/> xx f(<f>, h) = (j>-^-h. (2) 

Front solutions, <j> = <f>(x) where x = x — ct, of (2) [14] propagate in a preferred direction 
dictated by the minimization of T . The speed of a front connecting 4>+(h) at x = ~°° to 
4>-{h) at x = oo is given by 

c = fi (h) = a[U{<t>-)-U{<j> + )l (3) 

where a(h) = 1/ (p'(h) 2 dx is positive [11]. For negative h values, U{<p-) > U(4>+) and 
the front moves in the positive x direction (c > 0) so as to increase that part of the system 
having lower energy. When h is positive the front propagates toward negative x values 
(c < 0). Notice that (2) remains invariant under the transformation x — > —x. Thus, in 
addition to a front connecting <f>+(h) to 4>-(h) and propagating, say, at positive speed c, 
there exists a symmetric front connecting <p-(h) to <p+(K) and propagating at a negative 
speed — c. In the following we refer to these symmetric solutions as representing the same 
type of front solution. 

Imagine now that h is not constant but, instead, a second field h = h(x, t) coupled to 
4> = (f)(x,t). A variety of physical, chemical and biological systems fall in that category. 
Bistable optical systems [15], crystal growth [16], autocatalytic chemical reactions [2,8], and 
predator-prey systems [17], are a few examples. For the present purposes it is sufficient to 
consider the simplest case where h(x,t) is a diffusive field that responds linearly to changes 
in the order parameter <fi(x,t). More specifically we assume the form [8] 

h t = eg((f),h) + Sh xx g(<f>, h) = 4> - a - a x h, (4) 



-4- 

where e > is the ratio, t^/th, between the time scales associated with the two fields 
and h, S = D^/D^ is the ratio of diffusion constants, and the coefficient a\ is positive. 
The combined system (2) and (4) (denoted hereafter by (2+4)) is no longer variational. 
Yet, it resembles the original system (2) in having, for a proper choice of ciq and ai, three 
stationary homogeneous solutions of which two are stable. The stable solutions correspond 
to the intersection points of the nullcline g(c/),h) = with the branches 4> = 4>±(h) and 
are denoted here by ((f>±,h±) (see Figs. 1). Unlike the variational system, however, the 
two stable states can be connected by more than one type of front solutions when e is 
sufficiently small [10-13]. Preparing the system at the lower state (0_, h-) and perturbing 
it locally so as to induce a transition to the upper branch = (f) + (h), yields a front 
propagating to the right: (</>(x), h(x)) {4>±i h±) as X ~^ Too, c > 0. If, on the other 
hand, the initial state is the upper one, (4>+, h + ), a perturbation that induces a transition 
to the lower branch <j) = <f>-(h) yields a front connecting the same asymptotic states but 
propagating to the left: (</>(%), h(x)) — > (<f>±,h±) as x ~^ Too, c < 0. The two fronts 
are not related by the symmetry x — > — x and, therefore, represent two different types of 
front solutions. 

The multiplicity of front solutions and the symmetry x — > — x of (2) imply that along 
with a front that transforms the lower state to the upper state (<f> + ,h + ), there 

exists another front (hereafter "back" ) propagating in the same direction that transforms 
the upper state back to the lower one. A combination of the two may yield a persistent 
localized structure, provided there exists a mechanism which binds the back to the front. 

To study the emergence of such a structure we consider the small e regime, e < 1, and 
assume a nondiffusive h field, or 5 = (allowing diffusion of h will not affect the results 
qualitatively as long as S is not too large). We then distinguish between front and back 
regions where h barely changes, and outer regions where cj) can be eliminated adiabatically, 
4> = 4>±(h) [8-10]. Imagine now a front transforming the down state into the 

upper one (<f>+,h+) and propagating to the right. The front speed is determined by the 
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local value of h: c = fi(h-) > 0. A back that follows the front (so as to form a single 
up-state domain) will be affected by the field h(x) that develops behind the front. This 
field, as we will show below, provides the binding force the front exerts on the back. In 
order to find that field we insert <j> = <j>+(h) in (4) to obtain a closed equation for h, and 
use the boundary conditions h(x°f) = h- and h(x) — *• h + as x ~^ — oo, where x°f denotes 
the front position in a frame moving at speed c. Solving for h we find 

h( X ) = (h- - h+)e e "<x-*V +h+ X < X% (5) 

where h± = 2(±1 — ao)/(l + 2ai) and k = {a\ + l/2)/c. In deriving (5) we used the linear 
approximation 4>±{h) ~ ±1 — h/2 valid for small \h\. 

For an up-state domain to become a localized structure of fixed size, the back speed 
should be equal to the front speed c. We thus require = — c, where h® = h(xb) 1S the 

local value of h at the back. This relation determines h®. Note that h\ must be positive, 
whereas h-, the value of h at the front, is negative. Inserting h = h® and x = x° in (5) 
we find for the size of the localized structure 

According to (6), a localized structure of fixed size exists for h + > h® (for h + < h® the 
front speed is larger than the back speed and the up-state domain expands indefinitely). 
We will now show that this structure is also stable to translational perturbations. To this 
end we represent the back in the form 

(j>b(x, t) = 4> h (x - Xb) + ecj>l(x ~ Xb, et), 

where Xb = Xb + Xb(tt) is the actual back position, and write the level of h at the back 
as hb = h(xb) = h® + eh\{et). Using these forms in (2), assuming h is constant across the 
back, we find 

d 2 x <Pl + cd x <Pl + (1 - 30g 2 )^ = h\- e- 1 ^, (7) 
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where the dot over Xb denotes differentiation with respect to t. Solvability of (7) requires 
the right hand side of that equation to be orthogonal to (dcfy/dx) exp(c%). This leads to 

Xb = Pehl (8) 

where (3 is a positive constant. Using (5) to evaluate h b = h(xb) and consequently eh\ we 
find from (8) 

i b =/3(h+ -h b )(l-e e "*>). (9) 

The linearization of (9) about Xb = gives the equation Xb = —tKf3(h + — h^)xb- Since 
k > 0, we conclude that for e <C 1, a stable localized structure is formed whenever h + > h®. 

We consider now the other extreme, e ^> 1. Adiabatic elimination of h reduces (2+4) 
to the variational equation 

(j>t = (1 - a^)<P - <f) 3 + ao/cn + (f> xx (10) 

with h = 4>/ai — ao/ai, where we assumed that 5 <^ e. Equation (10) is equivalent to (2) 
and, consequently, has only one type of front solution connecting the two states (4>±, h±). 
When ao = the system (2+4) has an odd symmetry about (0, h) = (0, 0). The two states 
(0±j h±) are equally stable and the front that connects them is stationary. When 5 = 
it becomes an exact solution of (2+4) that exists for all e values. Fig. la shows a phase 
portrait of this front solution in the ((f), h) plane. It amounts to a straight diagonal line, 
h = a^V; connecting the two states. We recall that for e 1 equations (2+4) admit 
two types of propagating solutions. One may therefore expect to find a bifurcation from a 
single to multiple front solutions as e is decreased [12]. 

To study this bifurcation we consider the symmetric model (ao = 0) with a non- 
diffusive h field and write a propagating front solution, (f> = (f> p (x — ct), h = h p (x — ct), as 
power series in c 

4> p = 4> s + C(f)i + C 2 (p2 + ••• 

( n ) 

hp = h s + chi + c h 2 + 
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where (<j) s , h s = a^fVs) is the stationary front solution. Expanding e as well, e = 
eo+cei+c 2 e2 + ..., and using these expansions in (2+4) (with ao = S = 0) we find solvability 
conditions, one at each order, that determine the coefficients eo, ei, .... Carrying out this 
perturbation scheme to third order in c we find eo = a^ 2 , t\ = and €2 < 0. These 
results imply a supercritical pitchfork bifurcation occurring at e c = eo = a^ 2 . Near the 
bifurcation the front speed scales like c ~ (e c — e) 1 / 2 . Numerical studies on (2+4) confirm 
these results. A numerically computed bifurcation diagram for the symmetric case ao = 
is shown in Fig. 2a. 

The leading order corrections in (11) take the form 0i =0, h\ = 4>' s . Using these 
forms in (11) we see that the difference between the stationary and the counter propagating 
solutions close to the bifurcation point is that in the latter the field h is translated to the 
right or to the left by an amount proportional to c. This translation breaks the odd 
symmetry of the stationary solution and gives rise to phase portraits deviating from the 
diagonal, h = a|fV, as shown in Fig. lb. The speed c can be directly related to this 
deviation by using (11): 

/oo 
ip(j)' s dx if) = h p - a^cpp, (12) 

-00 

where a s = 1/ j^^i&'s) 2 dx 1S positive. Thus a front connecting (0 + , h + ) at x = — 00 to 
(<^>_, h-) at x = 00 propagates to the right (c > 0) when the deviation from the diagonal 
is negative (ifj < 0) and to the left when the deviation is positive (ifj > 0). As before, we 
refer to the latter case as describing a back. 

A significant distinction between the stationary front solution <j) = <f>s, h = a^fVs 
(in the symmetric model) and the two propagating solutions <j> = S , h = a^Vs ± c<p' s 
that bifurcate at e = e c can be made by looking at the phase 9 = arctan(/i/0). Across the 
stationary front the phase remains constant everywhere except for the core where it suffers 
a jump of 7r as the fields cf> and h vanish and change sign. Across a propagating front, on 
the other hand, the phase smoothly rotates by n keeping the modulus (0 2 + /i 2 ) 1 / 2 nonzero. 
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Similar types of walls have been found in the context of anisotropic ferromagnets [18] and, 
recently, in a nonvariational model describing a periodically forced oscillating medium [19] 
(see also Refs. [20,21]). They are referred to as Ising walls when the phase is singular and 
as Bloch walls when the phase rotates smoothly. 

The symmetric model is nongeneric unless the relevant physics dictates an odd sym- 
metry. Unfolding the symmetric case by allowing non zero values of ao yields a bifurcation 
diagram as shown in Fig. 26. In that case, front multiplicity arises by a saddle node bifurca- 
tion occurring at e sn (ao) < e c . For small ao we find the scaling behavior, e c — e sn ~ |a.o | 1//2 . 

In the small e regime we could show, using leading order singular perturbation anal- 
ysis, that front multiplicity gives rise to persistent patterns. We postpone the analogous 
analysis for higher e values to a subsequent study and present here, instead, results of 
direct numerical integration of (2+4) (with 6 = 0). Our observations are summerized in a 
phase diagram shown in Fig. 3 (see also Fig. 6 of Ref. [12]). To the right of the solid curve 
e = e sn (ao) only one type of front solution exists and, as we expect, initial patterns decay 
toward a uniform state. We note that along the line ao = the decay can be extremly 
slow and unnoticible in practical situations [22] . As we cross the solid curve by decreasing 
e below e sn (ao) front multiplicity arises. For some range of e initial patterns still decay 
toward a uniform state. Persistent patterns, in the form of stable solitary and periodic 
traveling waves, appear below a second critical value of e, that is, in the region to the left 
of the dashed curve. 

We have presented here a mechanism of pattern formation in systems undergoing a 
bifurcation from Ising to Bloch type fronts. A key ingredient in this mechanism is the 
coexistence of stable, counter propagating front solutions connecting the same asymptotic 
states as |x| — ► oo. Coexistence of such fronts cannot occur in variational systems having 
free energies. The mechanism can be tested in bistable chemical reactions [2] . 
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Figure Captions: 

Figure 1: Phase portraits of front solutions connecting the (0 + , h + ) state at x = — oo to 
the ((f)-, h-) state at x = 00 • The light colored curves are the nullclines / = and g = 
and the dark colored curves are the numerically computed trajectories. The computational 
parameters are: a. e = 1.0, 5 = 0, a\ = 2.0, ao = 0. b. e = .2, 5 = 0, ai = 2.0, ao = 0. 
c. 

Figure 2: Bifurcation diagrams of front solutions. The dots are data points representing 
the speed of the different types of stable front solutions that exist for each value of e. 
a. The symmetric case (ao = 0). b. The nonsymmetric case (ao = .1). 
Figure 3: Phase diagram in the e — a plane. For the region to the right of the solid curve 
only one type of front solution exists and initial patterns do not persist. In the region 
between the solid and the dashed curves multiple stable fronts coexist but patterns still 
decay toward a uniform state. For the region to the left of the dashed curve initial patterns 
evolve toward persistent patterns in the form of stable traveling waves. Computational 
parameters are : S = 0, a± = 2.0. 



ro 



c 

CD 

ro 



o 



CO -o 



